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DEPENDENCE OF ADSORPTION PROPERTIES ON SURFACE STRUCTURE 


FOR BODY-CENTERED-CUBIC SUBSTRATES 
by Robert J. Bacigalupi and Harold E. Neustadter 
Lewis Research Center 


SUMMARY 

The Lennard-Jones 6-12 atom interaction potential is applied to the cal- 
culation of adsorption energy of an atom onto as many as 132 sites on a unit 
cell area of each of the eight highest surface density planes of a body- 
centered-cubic substrate. From the calculations, topographical maps of adsorp- 
tion energy on a unit cell area of each plane for various adsorbate-adsorbent 
combinations are plotted. Normalized values of maximum adsorption energies and 
minimum surface diffusion activation energies for all cases are tabulated. The 
results are compared with available experimental data concerning adsorption of 
alkali metals, alkali earth metals, and inert gases on transition metal sub- 
strates and are in good agreement. Application of the results to transition- 
metal - transition-me cal combination predicts that the 110 surface is the 
lowest energy configuration for a body-centered-cubic crystal. 


INTRODUCTION 

Adsorption on an atomic scale gained recognition with the advent of the 
field-ion microscope and field-emission microscope, and it gained in importance 
proportionally with the fields of ultrahigh vacuum, gas-filled thermionic con- 
verters and the contact-ionization engine. Through the use of field-ion and 
field-emission microscopy, it has been observed that adsorption properties are 
strongly dependent on the atomic arrangement of the substrate (in some cases 
more so than on the materials involved) . Forces governing adsorption of many 
species on metals or semiconductors have been found to be of the form of dis- 
persion forces such as van der Waal forces (ref. l) . By considering dispersion 
force interactions of an adsorbing atom with the nearest neighbors of the sub- 
strate, one can verify that the substrate structure grossly affects the adsorp- 
tion properties (refs. 2 and 3). It is the intent of this report to calculate 
the interaction potential encountered by a single atom at any point on various 
ideal single crystal planes of the body-centered-cubic structure. Furthermore, 
it is intended to deduce from such potential surfaces some general properties 
of physical adsorption and specific properties, such as adsorption energy and 
diffusion activation energy, and to compare these results with experimental ob- 
servations . 


The Lennard- Jones 6-12 potential is used as a model. Its applicability 
to adsorption of a nonpolar^- atom on a metal is discussed in the following sec- 
tion. The general mathematical techniques employed are outlined in the section 
PROCEDURE, while details concerning the performance of the specific lattice sum- 
mations for the body-centered-cubic structure are contained in appendix A. The 
last section compares the results with the scant experimental data and discusses 
the application of the results to known adsorption characteristics of practical 
surfaces . 


MODEL 

The Lennard- Jones potential is (ref. 2) 

E = 4E te 6+ (?f] (1) 

where E is the atom-atom interaction energy, e is the depth of the energy 
well, a is that finite value of r for which E is zero, and r is the in- 
ternuclear distance. The attractive portion of equation (l) is proportional to 
the inverse sixth power and represents a long-range van der Waals interaction, 
while the repulsive contribution is approximated by an inverse twelfth power de- 
pendence. 


Application of this atom-atom interaction potential to adsorption on a 
solid assumes additivity of the Interaction of the adsorbate with each atom of 
the adsorbent or 

^ 4E Zte 6+ (jf] < 2 > 

i 

where 9 is the atom-metal interaction energy and r^ is the distance from 
the adsorbed atom the i^* 1 atom of the substrate. Thus, effects arising specif- 
ically from the structure of the substrate are contained in the lattice summa- 
tions, X) r ~ 6 and r i"^ whereas effects arising from the atomic properties 
i i 

of the constituent atoms are contained in the parameters a and e. Appen- 
dix B contains a discussion of the physical meanings of both a and e. 

Equation (2) neglects any contribution to the interaction by free electrons 
in the solid; however, a study by Pierotti and Halsey (ref. 4) on the inter- 
action of krypton atoms with metals has shown that theories which include the 
free electron contribution of the metal do not give significantly improved re- 
sults as opposed to theories in which the free electron contributions are ne- 


-*-All atoms have instantaneous dipole moments that are the primary source 
for the long range interactions between atoms. An atom is considered nonpolar 
if when isolated its net dipole moment is zero over a period of time. 
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glected. Neustadter, Luke, and Sheahan (ref. 5) have further shown that very 
low coverage adsorption of various alkali metal atoms on tungsten can be satis- 
factorily explained by using equation (2). Of course, the validity of the pro- 
posed model can best be tested by comparing the calculated results with experi- 
mental data. 


PROCEDURE 


In order to preserve the general applicability of this work to adsorption 
on any substrate, equation (2) is rewritten as 


JL 

4e 



( 3 ) 


where d^ is r^/a Q , and a Q is the lattice constant of the substrate. Equa- 
tion (3) is treated as follows: For each surface orientation considered, a co- 

ordinate system is chosen such that the x,y plane passes through the centers 
of all the surface atoms and z is positive along the outward normal. Clearly 
qj/4e is a function of x, y, and z. 


As shown in sketch (c) of appendix A (p. 9) each unit cell is broken up 
into a grid of as many as 132 prospective adsorption sites where the number of 
sites is chosen so that each unit cell will have approximately the same density 


of prospective sites. At any one of these sites (x i ,y i ) the sums 




and 


F.AT12 are evaluated as functions 



(a) 



(b) 


of z. Using these sums in equa- 
tions (3) and plotting cp/4e as a 
function of z yield, for example, 
sketch (a) (which will depend on sur- 
face orientation, site, and cr/a 0 
for its specific form). 

The assertion is made that, if 
adsorption were to occur over site 
x i ,y i , the nucleus of the particle 

would be located at z m p n , and the 
bond energy would be cp m in/ 4e (where 
the quantities are as defined in 
sketch (a)). The derivative of cp/4e 
is taken (numerically) with respect 
to z , and Tmin/ an A z min are 
calculated for a specific a/a 0 (see 
appendix A). This entire process is 
now repeated for each predetermined 
location xj,yj on a unit cell of 
the adsorbent. Then cp m in/ is 
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plotted as a function of surface coordinates x and y. An example of such a 
plot is shown in sketch (h) for the 110 plane and a/ a 0 = 1.00 (corresponding 
to barium on tungsten). Here one unit cell of surface area is represented with 
the locations indicated for the substrate atoms nearest the surface. Along 
each of five equally spaced lines of constant y,cp m ^ n /4e is plotted against x 
with the ordinate cp^ j-n/ 4e in the y-direction. 

The minimum value of cp m i n /4e is read graphically and is listed in 
table I. These values are the normalized heat of adsorption 3>/4e for the pre- 
ferred adsorption site, for a specific orientation, and for a specific value 
of a/a Q . The quantity 4e is also determined graphically by identifying 
it with the least amount of energy that must be supplied to an adsorbed atom to 
enable it to move from one preferred adsorption site to another, where the 
assumed path was not necessarily a straight line. In the 110, 210, and 310 
cases, the topographical potential plots were obtained by manually cross- 
plotting the data for cp m in/ 4e against position. In the remaining cases, the 
topographical potential plots were generated directly by a Fourier interpolation 
7090 program from the original digital output. This procedure was made possi- 
ble by using the fact that the potential is periodic along any line parallel to 
a border of the unit surface cell. This entire procedure is repeated for a 
range of values of a/a 0 and the remaining substrate orientations. For the 
specific summation technique applied to the body-centered-cubic structure see 
appendix A. 


DISCUSSION OF RESULTS 

The minimum values of cp/4e and E^/ie are listed in table I. The orig- 
inal plots of the potential variation, such as illustrated in sketch (b) , were 
replotted as topographical maps showing lines of equal adsorption energy . The 
eight highest surface density planes were considered (110, 100, 211, 310, 111, 
321, 411, and 210 in order of decreasing density) and are presented in fig- 
ures 1 to 8, respectively. Included in the topographical maps are the locations 
of the substrate atoms in each unit cell. The most probable migration direction 
is assumed to be path of least resistance (minimum E^/ie) connecting two equiv- 
alent adsorption sites. 

The values of o/a Q used in the calculations were 0.80, 0.894, 1.00, 

1.123, and 1.20. Table II lists the values of cr/a Q corresponding to various 
adsorbent -adsorbate combinations. The most meaningful interpretation of these 
data results from considering the comparative effects of substrate surface 
structure (crystal orientation) on adsorption and diffusion for a given combi- 
nation of adsorbent and adsorbate (i.e., for a given a/a Q ). 

The Lennard-Jones potential model corresponds most closely to inert gas 
atoms adsorbing on a semiconductor or an insulator when the coverage is low 
enough to eliminate any interaction of adsorbate atoms with one another. This 
limit of coverage is shown experimentally in some cases to be as high as 0.7 
(e.g., thorium, ref. 6), thus allowing the calculations to be useful in 
practical systems at coverages other than zero. The Lennard-Jones potential 
model has been used in calculating adsorption energies of various metals on 
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tungsten (ref. 7) and, in particular, adsorption energies of alkali metals on 
tungsten (ref. 5) . Agreement of the results with experimental data implies 
that this method may also be applied to other metal -metal systems. Substrates 
that are compared are assumed to be atomically smooth ideal crystal planes with 
no dislocations, steps, impurities, missing atoms, or surface relaxations. The 
experimental situation most closely approaching these assumptions exists in the 
field-emission microscope in which the tip is an annealed single crystal with a 
minimum of imperfections and on which there are small areas of atomically smooth 
planes . 

In emission-microscope experiments by Utsugi and Gomer (ref. 8) and 
Drechsler (ref. 9) who studied the adsorption of barium on tungsten, the acti- 
vation energies for diffusion E d on the 110 and 321 planes were measured and 
were in good agreement with the calculated results as follows: 

Measured: 


(ll°)E d Q.4! 

(321 )E d 0.83 


0.49 (Ref. 8) 


( 110 )E d 0>20 

("321)E d ~ 0.65 


0.31 (Ref. 9) 


Calculated: 

(H0)E d /46 ^ Q 155 _ Q 39 
(32l)E d /4c 0.345 

Experimental data such as these may be used to evaluate e empirically, which 
in turn may be used to obtain quantitative values for and 0 for each 

plane of a given system (i.e., given a/a Q ) . An alternative theoretical expres- 
sion for e is presented in appendix B. 

Calculated adsorption energies from table I (with a/a Q = 1.00 for xenon 
on tungsten and molybdenum) agree favorably with emission microscope experi- 
ments by Ehrlich and Hudda (ref. 10) , who studied the adsorption of xenon on 
tungsten and molybdenum. They found stability of the adsorbate highest on the 
210, 310, and 611 planes, lower on the 211 plane, and lowest on the 100 and 110 
planes at temperatures where diffusion is not important. This remarkable agree- 
ment is less surprising when one considers that the model is best suited to 
inert gas interactions. 

Adsorption of cesium on transition metals is of considerable interest in 
those applications related to electron and ion emission. Unfortunately, the 
adsorption energy of cesium on various planes of metals is not as well docu- 
mented as that of inert gases on metals. It has been repeatedly observed (e.g., 
see ref. ll) that the 110 plane of both molybdenum and tungsten has a higher 
current density than any other plane at a given substrate temperature and ce- 
sium arrival rate (in the region of low to intermediate coverages). Since this 
plane has the highest vacuum work function, these observations have been inter- 
preted as implying a maximum stability of cesium on the 110 surface, whereas 
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the calculations predict that cesium should be least stable on the 110 plane of 
tungsten ( a/a Q = 1.123). However, one may attribute increased emission density 
from a partly covered plane to any one or combination of the following three 
factors : 

(1) High stability of the adsorbate on the plane in question, which de- 

pends exponentially on 0 and E^ 

(2) High number of sites available for adsorption 

(3) Large effective dipole moment of an adsorbed particle 

The 110 plane has the highest density of sites available of all body- 
centered-cubic planes. Good and Muller (ref. 6) have shown that for a given 
surface density of adatoms the 110 showed by far the greatest decrease in work 
function of all body-centered-cubic planes. The relatively large emission den- 
sity from the 110 plane of tungsten partly covered with cesium is not suffi- 
ciently definitive of adsorption energies to provide an adequate test of the 
Lennard- Jones potential model. 

In general, the calculations substantiate the prevalent model of surface 
adsorption in which discrete adsorption sites exist and in which diffusion is 
defined as the motion of atoms between well-defined sites. Also the calcula- 
tions verify the expectation that these sites exist at wells in the surface 
where the outer electrons of adsorbate participate in binding with a maximum 
number of substrate atoms. 

Additional verification of the validity of the calculations is realized in 
examining self-adsorption and self -diffusion. The lowest energy surface config- 
uration of a body-centered-cubic crystal of tungsten (refs. 12 and 13) was shown 
experimentally to be the 110 surface, where the experimental situation closely 
approaches that assumed in the calculations. Recognizing the relatively low 
values of ®/4e and E^/le a_ t cr/a Q = 0.80 (cr/a 0 = 0.786 for self-adsorption) 
for the 110, 100, and 211 planes (in that order) one can see that the calcula- 
tions predict that, for a body-centered-cubic crystal annealed in vacuum, the 
final surface would consist of 110 planes, while, in the transition, the 100 
and 211 planes have a high probability of existing. It must be noted here that 
even though the physical mechanism of metallic binding may not be described ex- 
actly by the Lennard-Jones potential, the structure-dependent adsorption prop- 
erties of such a system can be readily predicted by using the Lennard-Jones po- 
tential summation. 


SUMMARY OF RESULTS 

The following four factors serve to support the applicability of the results 
to practical adsorption phenomenon: 

1. The general picture of adsorption, as observed experimentally, is ver- 
ified because atomically rough surfaces tend to have higher heats of adsorption 
than closer packed surfaces. 
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2. More specif ically, the order of heats of adsorption for xenon on var- 
ious planes of tungsten is the same as that observed in field-emission- 
microscope experiments. 

3. The equilibrium surface of body-centered cubic crystals is the 110 
plane determined both by experiment and calculation. 

4. The relative diffusion activation energy for barium on tungsten is in 
good agreement with measured values. The possible extension of these calcula- 
tions to the adsorption properties of other structures (e.g. , face-centered- 
cubic, diamond, etc.) and of single-crystal surfaces with various imperfections 
(e.g., missing atoms and steps) can readily be seen. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, August 6, 1965. 
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APPENDIX A 


DETAILS OF THE SUMMATION FOR BODY-CENTERED-CUBIC SUBSTRATES 

The location of all atoms in a body-centered-cubic crystal with respect to 
the coordinate system described previously has been recorded by Bacigalupi 
(ref. 14) for twenty orientations. This_ t£cbnique can be adapted to the present 
problem to define the four quantities S, X, V, and d, where S locates the 
nearest surface atom to any surface atom, X locates the nearest surface neigh- 
bors in the x -direction to any surface 
atom, V connects any surface atom to 
the closest atom in the next-to-surface 
plane, and d is the interplanar dis- 
tance measured along the z-axis. The 
vectors X and S form a unit cell 
of the surface as shown in sketch (c). 

With these definitions and equa- 
tion ( Al) all the atoms in the crystal 
with respect to the surface plane can 
be located 

x i * X o * A|V(X)| + n|X| + C 

Y i = Y o + Y i|T( Y )| + 1|S(Y)| 

C - [y, - Y 0 - i^ODi ] j§[|j 

where X^,Y^ is the location of any atoms referred to the surface plane (x,y 
plane), X Q ,Y 0 is the starting point and location of the upper left corner of 
the unit cell, 1^ is the plane in which the i^* 1 atom exists (0 for surface 
plane; 1 for next-to-surface plane, etc.) (x) or (Y) indicates the component of 
the particular vector, and n and m are integers. 

In performing the summations, it was necessary to establish a cutoff for 
the computer. Using an integral approximation to determine the contribution of 
those terms omitted from the summation leads to the result that accuracy to 
within 1 percent is achieved when the summation is extended over all the atoms 
within seven a Q of the site being considered. Therefore equation (A2) is in- 
troduced to establish the cutoff: 




X Q + 7.0 + X| > X j _ > X Q - 7.0 


Y o +7.0>Y 1 >Y o -7.0+ |S(Y)J 


The prospective adsorption sites for which the summations were performed 
are located at the intersections of five to ten equidistant lines parallel 
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to X and eleven to twenty equidistant lines parallel to S, as shown in 
sketch (c). The distance rj_ from an adsorbing atom to the ith atom in the 
crystal can now be evaluated as a function of z. At each site, d 1 (defined 
by eq. (A3)) 

a i - [( X 1 - x o) a + (*i - *o) 2 + (M + z ) 2 ] l/2 (A3) 

was evaluated for each of sixty evenly spaced values of z ranging from 0.05 
a Q to 2.1 a . The corresponding values of dj_ were used in the summation in 
equation (3). The minimum value of cp/4e obtained from equation (3) for a 
given site was used along with one value of cp/4e in the negative z-direction 
and three values in the positive z-direction. A fourth-order equation was 
fitted to these five points , and a min imum was determined. The minimum was 
further corrected by using the Newton-Raphson method. Sketch (b) is an ex- 
ample of these minima plotted as a function of their position for the 110 
plane for cr/a 0 = 1.00. 
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APPENDIX B 


EVALUATION OF a/a Q AND e 

The Lennard-Jones potential (ref. 2) for interactions between two atoms is 


E = 4e 


[-©••iff] 


( 1 ) 



which graphically is of the form 
shown in sketch (d), where a is 
that finite value of r for which 
E = 0, e is the interaction bond 
energy of the two atoms, and r Q 
is their internuclear equilibrium 
distance . 


To evaluate a note that 
is that value of r for which E 
is a minimum. Therefore, a can 
be related to r Q by setting 
(dE/dr) r = 0. This gives 


o 


= 0. T1 

a = r 0 2"l/6. Since r can be ap- 
proximated as the sum of the hard- 
sphere radii of the two interacting 


atoms, values can now be calculated for c. Table II lists cr/a 0 for various 
combinations of adsorbates and adsorbents with body-centered-cubic structures. 


Various theories have been developed to evaluate explicitly the attractive 
long-range interaction energy for two dissimilar atoms. The best results seem 
to be obtained (ref. 4) from the theory developed by Kirkwood (ref. 15) and 
Muller (ref. 16) . They obtained 


E (attractive) 


6mc 2 a A a M 
r 6 a* 

X A 


(Bl) 


where a is the electronic polarizability, X is the diamagnetic susceptibility, 
m is the electron rest mass, c is the speed of light in vacuum, the subscript 
A refers to the adsorbed atom, and the subscript M refers to the substrate 

metal atoms . 


Comparison of equation (Bl) with the attractive term in the Lennard-Jones 
potential allows evaluation of e from 


4ea 6 = 


6mc^ 


a A a M 


^A + ^ 

X A % 


(B2) 
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TABLE I. - NORMALIZED VALUES OF HEATS OF ADSORPTION AND DIFFUSION ACTIVATION 
ENERGY FOR VARIOUS ADSORBATES ON BODY -CENTERED -CUBIC SUBSTRATES 


Plane 

(body 

centered 

cubic) 

Q 

O 

0.80 

0.894 

1.00 

1.123 

1.20 

$/4e 

%/« 

$/4e 

E d /« 

$/4e 

Ea/4e 

$/4e 

V 4e 

$/4e 

V« 

110 

1.29 

0.163 

1.54 

0.152 

1.90 

0.135 

2.41 

0.110 

2.79 

0.093 

100 

1.581 

0.586 

1.732 

0.514 

1.991 

0.449 

2.401 

0.382 

2.719 

0.343 

211 

1.761 

0.308 

2.045 

0.387 

2.322 

0.346 

2.707' 

0.262 

3.009 

0.216 

210 

1.84 

0.45 

2.02 

0.28 

2.39 

0.18 

2.97 

0.12 

3.35 

0.10 

111 


0.680 

1.933 

0.480 

2.125 

0.397 

2.484 

0.340 

2.772 


310 


0.56 

2.10 

0.41 

2.34 

0.31 

2.78 

0.17 

3.04 

0.07 

321 

1.815 

0.310 

2.162 

0.422 

2.509 

0.345 

3.025 

0.275 

3.417 

0.231 

411 

1.930 

0.920 

2.221 

0.925 

2.466 

0.783 

2.802 

0.640 

3.071 

0.510 
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TABLE II. - VALUES OF a/ a Q FOR VARIOUS ADSORBATES 
ON BODY -CENTERED CUBIC SUBSTRATES 3, 


Adsorbate 

Substrate 


Vanadium 

Chromium 

Iron 

Tantalum 

Niobium 

Molybdenum 

Tungsten 

Neon 

0.864 

0.890 

0.891 

0.825 

0.828 

0.846 

0.844 

Argon 

0.955 

0.985 

0.987 

0.909 

0.911 

0.934 

0.932 

Krypton 

0.981 

1.013 

1.015 

0.933 

0.935 

0.960 

0.957 

Xenon 

1.040 

1.075 

1.077 

0.987 

0.990 

1.017 

1.013 

Lithium 

0.849 

0.874 

0.875 

0.812 

0.813 

0.832 

0.830 

Sodium 

0.951 

0.981 

0.982 

0.905 

0.908 

0.931 

0.927 

Potassium 

1.084 

1.122 

1.124 

1.027 

1.029 

1.064 

1.061 

Rubidium 

1.227 

1.162 

1.164 

1.063 

1.066 

1.096 

1.092 

Cesium 

1.178 

1.220 

1.224 

1.114 

1.117 

1.150 

1.123 

Thorium 

0.921 

0.950 

0.951 

0.878 

0.880 

0.902 

0.894 

Barium 

1.045 

1.080 

1.083 

0.992 

0.994 

1.021 

1.018 

Strontium 

1.017 

1.050 

1.053 

0.966 

0.968 

0.994 

0.991 

Self- 

adsorption 

0.784 

0.785 

0.784 

0.786 

0.787 

0.787 

0.785 


a Based on atomic radii obtained from ref. 17, p. 403. 
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Figure 1. - Topographical map of normalized adsorption energy on 110 plane. 
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Figure 2. - Topographical map of normal- 
ized adsorption energy on 100 plane. 
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Figure 5. - Topographical map of normalized adsorption energy on 111 plane. 
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(c) o la Q = 1. 20. 

Figure 6. - Topographical map of normalized adsorption energy on 321 plane. 







(c! a/a 0 > 1. 20. 

Figure 7. - Topographical map of normalized adsorption energy on 411 plane 
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"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 
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